Spin and Conductance-Peak-Spacing Distributions in Large Quantum Dots: 
A Density Functional Theory Study 
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We use spin-density-functional theory to study the spacing between conductance peaks and the ground-state 
spin of 2D model quantum dots with up to 200 electrons. Distributions for different ranges of electron number 
are obtained in both symmetric and asymmetric potentials. The even/odd effect is pronounced for small sym- 
metric dots but vanishes for large asymmetric ones, suggesting substantially stronger interaction effects than 
expected. The fraction of high-spin ground states is remarkably large. 
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The interplay of quantum mechanical interference and 
electron-electron interactions is a current theme in many areas 
of solid-state physics: the 2D metal-insulator transition, inter- 
action corrections in mesoscopic systems, and efforts toward 
solid-state quantum computing, for instance. A semiconduc- 
tor quantum dot (QD) [fjj [p ~ a nano-device in which electron 
motion is quantized in all three dimensions - is a particularly 
simple system in which to study this interplay. In Coulomb 
blockade experiments in the electron tunneling regime, the 
conductance through the dot varies strongly as a function of 
gate voltage, forming a series of sharp peaks. For closed dots 
at low temperature, both the positions and heights of the peaks 
encode information about the dot's ground state. In particu- 
lar, the spacing between adjacent conductance peaks is pro- 
portional to the second difference of the ground state energy 
with respect to electron number N, A 2 E(N)=E gs (N + 1)+ 
E gs (N - 1) - 2E gs (N), which is often called the addition 
energy. Furthermore, the ground state spin of the QD can be 
inferred from the shift in position of the conductance peaks 
upon applying a magnetic field. 

The addition energy varies because of changing interfer- 
ence conditions either as N changes or from dot to dot, lead- 
ing to a conductance-peak-spacing distribution. Previous the- 
oretical work addressing this distribution can be divided into 
roughly two types: First, computational approaches addressed 
small dots with randomly disordered potentials - both ex- 
act diagonalization J3, EL 5|] and self-consistent field meth- 
ods (Hartree-Fock ^T^f^^, [[(J or density functional the- 
ory mU, n~2|). Second, a semi-analytic treatment of large 



dots was developed based on general statistical assumptions 
[§, [l3[ |l4t |l5, 16, ivt|: An important contribution to the varia- 
tion comes from the single-particle energy; to treat this for ir- 
regular quantum dots, one assumes that the single-particle dy- 
namics is classically chaotic, and so the single-particle quan- 
tum properties can be described by random matrix theory 
(RMT) [|J |l8l ], A random-phase approximation (RPA) treat- 
ment of the screened electron-electron interaction was then 
combined with such an RMT description of single-particle 
states to describe dots with large N. 

The results of these two approaches are quite different. 



First, for the zero temperature peak-spacing, the small dot cal- 
culations yield Gaussian-like distributions while the large N 
results are non-Gaussian. Second, spin degeneracy causes a 
significant "even/odd effect" in the large dot approach: the 
distribution for N even is very different from that for N odd. 
Third, with regard to the ground state spin [|l2 , 19 2^, [H] |22|] , 
the small N calculations find an enhancement of the low spin 
states compared to the large QD approach. 

Experimental work to date ^ [23] ^4], [Z5] ^(J has unfortu- 
nately failed to probe the ground state addition energy distri- 
bution or ground state spin of generic systems. In the more 
recent experiments, either temperature obscured the ground 
state properties [|l5 , 16 , 24], |6|] or the dot was regular in shape 



Our aim here is to bridge the gap between the two theo- 
retical approaches and in so doing highlight the need for more 
experiments. We have used the Kohn-Sham (KS) spin density- 
functional theory (SDFT) [|7| to study both the peak-spacing 
and the spin distribution for 2D model QDs. Using an efficient 
algorithm [[35]], we obtain statistics with N up to 200. This is 
the first calculation, as far as we know, for large realistic QDs. 

Our primary result is that the effective electron-electron in- 
teraction that emerges is substantially stronger than that pre- 
dicted from the RPA-RMT treatment. The evidence is two- 
fold: (1) the addition energy distribution is Gaussian-like with 
no discernible even/odd effect, and (2) the probability of hav- 
ing a high spin state (S > 1) is larger than the maximum 
possible from RPA-RMT. 

In SDFT, the energy of the interacting system is expressed 
as a functional of the spin up (a) and down (/?) electron den- 
sities, p a (r) and pp(r). (Effective atomic units are used: 
for GaAs-AlGaAs QDs, the values are 10.08 meV for en- 
ergy and 10.95 nm for length.) The functional is built out 
of four parts. Two are trivial - the response to the external 
potential V ex t and the Coulomb energy. The kinetic energy 
of a non-interacting reference system is included explicitly: 

T s [p a , P p] = E^Wwl - ^v 2 |v^) with Pa = £v |^| 2 

and a = a, f3. The final term is, of course, the exchange- 
correlation energy E xc [p a ,pp], Thus the KS-SDFT func- 
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FIG. 1 : Addition energy as a function of electron number, (a) Raw 
addition energy data (solid) and its polynomial fit (dashed) to remove 
the change in classical charging energy, (b) Addition energies scaled 
by the mean level-spacing after removing smooth part. (A = 0.6 and 

7 = 0.) 



tional is 



E[p a , pp] = T s [p a , pp] + J V ext {v)p{v)dv 
P \ r)p{r \ vdv' +E xc [p a ,pp]. 



(1) 



The ground state energy is obtained by minimizing the func- 
tional with respect to the spin densities under the constraints 
J p a (r)dr = N a . We use the standard local spin density 
approximation (LSDA) for E xc which works well in various 
semiconductors; in particular, we use Tanatar and Ceperley's 
2D parametrization p8|]. For 2D clean QDs, comparisons 
with quantum Monte-Carlo calculations for N < 8 and in- 
teraction strength parameter r s < 8.0 have shown that LSDA 
works well for both the ground state spin and energy 
We use a quartic potential to model 2D QD systems, 



x 4 

V ext [r) = a[— + bi/ - 2AxV + l(x V ~ xy 2 )i 



(2) 



Both the classical dynamics and the single-particle quantum 
mechanics at 7 = have been studied in detail [|T]]: the sys- 
tem evolves continuously from integrable to fully chaotic as 
A changes from to 1 . The parameter 7 breaks the four-fold 
symmetry. The prefactor is a = 1.0 x 10 -4 ; this allows the 
electrons to spread so that the interaction strength parameter, 
r s , is about 1.5, close to experimental conditions. 

For a given V ex t, we calculate the total energy at several val- 
ues of electron number N and total spin S. Selecting the min- 
imum energy determines the ground state energy E gs and spin 
S gs as a function of N. The addition energy is then calculated 
as the second difference of E gs {N); an example is shown in 
Fig. ]j](a). To obtain good statistics, in both the symmetric and 
asymmetric cases we calculate five sets of data with different 
parameters [[36]). A correlation analysis shows that both the 
single-particle level spacing (SPLS) and addition energy from 
the different sets are statistically independent. 

Since the slow decrease in A2E(N) is a classical effect - 
the increasing capacitance as the dot becomes bigger - we 




FIG. 2: Distributions of peak-spacing, s, for even (solid) and odd 
(dashed) N. The columns correspond to the different ranges of elec- 
tron number indicated, while the first row [(a)-(c)] is for a symmetric 
external potential and the second for asymmetric [(d)-(f)]. The inset 
of each figure shows the corresponding spin distribution for N even 
(sold) and odd (skeleton). Both the lack of even/odd effect and the 
large spin in panel (f) is striking and indicates an unexpectedly large 
effective electron-electron interaction. (A sliding window is used 
in estimating the probability density, yielding a smooth curve rather 
than a histogram; each curve is made from 150 data points using a 
Gaussian window of width 0.3.) 



remove it by fitting a polynomial to find the smooth part 
(A2E(N)). To compare with experiments, the addition en- 
ergy is scaled by the mean-level spacing found from the av- 
erage electron density n through A = 2nh 2 (n) / 'm* N . The 
resulting dimensionless spacing is denoted s = [A2E(N) — 
{A 2 E(N))]/ A; see Fig. [jjb). Note that the typical scale of s 
is 1 ; that is, fluctuations of the addition energy are on the scale 
of the single-particle mean-level spacing. 

We find the distribution of s for even and odd N in three 
ranges of electron number, A^ = 20-80, 80-140, and 140-200, 
in both symmetric and asymmetric external potentials. (For 
the smallest, mean, and largest values of N, r s is 2.5, 1.6, and 
1 .4, respectively.) The main features of the results, shown in 
Fig. |] are as follows. (1) Even/odd: There is a difference in 
the distribution for N even or odd in the small A" range, and 
the parity effect in the symmetric case is more pronounced 
than for asymmetric potentials. But there is a striking absence 
of even/odd effect in the asymmetric large A^ case [panel (f)]. 
(2) Shape: The distributions are Gaussian-like. This extends 
the previous small A^ results, disagrees with the RPA-RMT re- 
sults for large AT, and is consistent with experiments [^4] ^6|] 
(where, however, the ground state properties are obscured by 
temperature effects JL5] , |T6|]). (3) Small/large: The spacing 
distributions for small A^ are different from those for large 
N. Since the experimental QDs generally involve tens to 
hundreds of electrons, one must be cautious in generalizing 
to large dots conclusions drawn from studying small dots. 
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(4) Symmetric/Asymmetric: Both the variance of the peak- 
spacing and the magnitude of the even/odd effect is larger in 
the symmetric case. 

The insets in Fig. || show histograms of the spin distribu- 
tion for the corresponding electron number ranges. A remark- 
able feature is the significantly higher fraction of high spin 
ground states that we find at large N than in either previous 
SDFT investigations of small disordered dots fll2| ] (for small 
N, we agree, of course, with previous results) or other in- 
vestigations. Especially in the asymmetric case, P(S = 1) is 
even higher than P(S = 0) for N — 80-200! Again we see 
a clear dependence on the electron number range: while the 
spin distribution for odd N changes little as N increases, the 
spin distribution for even N is quite different for small and 
large N, Comparing spin distributions in the presence and ab- 
sence of symmetry, we see a larger high-spin fraction in the 
asymmetric case for all three ranges of N. 

This last result is at first surprising: generally one expects 
increased symmetry to increase the spin - as in Hund's rule 
for the spin of atoms. We can use random matrix theory, how- 
ever, to show that there are competing effects here - namely 
the statistics of the eigenenergies vs. the statistics of the eigen- 
functions. For simplicity we consider the simplest RMT: a 
two-electron two-level model in which we neglect spatial cor- 
relations beyond a wavelength. In a Hartree-Fock framework 
and under the assumption that the orbitals remain the same for 
different spin configurations, the energy difference between 
the singlet (S = 0) and triplet (S = 1) states is 



SE = E s =i - E s=0 = Se + (J 12 - K ls 



Jn 



(3) 



where Se is the single-particle level spacing, Jy = 
J drdr'\ijji(r)\ 2 v scr (r, r')|-0j(r)| 2 are the Coulomb energies, 
andiCi2 = J drtfrV5'(r)V;|(r> acr (r,r')V'i(r')V'2(r) is the 
exchange energy. Here instead of using the bare Coulomb 
interaction, we use the screened potential v scr (r, r') in or- 
der to implicitly account for the other electrons. Qualita- 
tively, the screened interaction is approximately zero-range, 
v scr (r, r') — » (AA/2)S(r — r') where A is the mean single- 
particle level spacing and A is the area. In this limit, the 
second term in Eq. (Q) vanishes, and the third term is pro- 
portional to the "inverse participation ratio" (IPR) defined by 
I = A J G?r|i/>(r)| 4 . Therefore, in the zero-range limit and con- 
sidering the time-reversibility of the system [B7fl, we have 



SE. 



zero — r;in<;< - 



= 5e- I A/3 . 



(4) 



There is clearly a competition here between the level spacing 
- a large spacing tends to decrease the spin - and the statistics 
of the wave functions - increased localization increases / and 
leads to a larger spin. 

To see the competition explicitly, we calculate the distribu- 
tion of the SPLS 5e and that of the IPR from top-level KS 
orbital energies and wave functions in both the symmetric 
and asymmetric cases. Results are shown in Fig. 0. While 
the symmetric case has a higher probability of small level- 
spacing, the mean value of the IPR is also smaller in the sym- 
metric case. Our overall result for the spin distribution - the 




FIG. 3: Distributions of the single-particle level spacing (scaled by 
the mean) and the IPR from Kohn-Sham calculations and from ran- 
dom matrix theory, (a) The distribution of SPLS calculated from top 
KS orbital energies for symmetric (solid) and asymmetric (dashed) 
case respectively, (b) The distribution of IPR calculated from top- 
level KS orbital wave functions for symmetric (solid) and asymmet- 
ric (dashed) case respectively, (c) The distributions of SPLS for a sin- 
gle GOE (dashed) and the superposition of four GOE (solid) accord- 
ing to RMT. (d) The distributions of IPR calculated from the eigen- 
vectors of real symmetric random matrix with dimension M — 20 
(solid) and M = 80 (dashed). 



decrease in probability of S = 1 upon introducing symmetry 
[Fig. 2(f)] - shows that the effect of wave function statistics 
is stronger in our system. 

The trend here is captured by the random matrix theory. For 
an asymmetric chaotic external potential, the distribution of 
the SPLS is that of the Gaussian orthogonal ensemble (GOE). 
In the symmetric chaotic case, however, the SPLS statistics is 
the superposition of four GOE's [^8[], one for each symmetry 
class. Fig. |3|(c) shows these two distributions. Clearly, the su- 
perposition greatly reduces the nearest-neighbor level repul- 
sion, which implies that spatial symmetry favors a high-spin 
ground state, in accordance with Hund's rule. 

For the wave function statistics, RMT suggests that the 
single-particle wave functions of classically chaotic systems 
are described by A/-dimensional random unit vectors [g, |32] ]. 
The shape of the resulting distribution of IPR depends on M, 
as shown in Fig. ||(d) for M — 20 and 80. These values of M 
are chosen to correspond to the number of independent orbital 
levels: approximately N/2 in the asymmetric case because of 
spin, and NJ 8 for the four-fold symmetric external potential. 
Symmetry reduces the effective M and so reduces the IPR, 
acting against a high-spin ground state. Hence there is a com- 
petition between the SPLS and the IPR. 

Note that the dependence of the IPR on M also explains 
the change in the spin distribution as electron number varies: 
The wave functions for small N have a smaller effective M, 
hence a smaller IPR. Thus a low-spin ground state is more fa- 
vorable. Similarly, a three-level analysis provides a qualitative 
explanation of odd N results. 

The results in Fig. || show, of course, that the calculated 
distributions of both the SPLS and the IPR agree with RMT 
only qualitatively. The fluctuation of the IPR is particularly 
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striking: it is much smaller for the KS wave functions than in 
RMT. We believe this is due to the neglect of spatial correla- 
tions in our very simple RMT. 

In summary, by studying a model 2D quantum dot with 
up to N = 200 electrons, we have found new phenomena. 
Both the statistics of ground state spin and the spacing be- 
tween conductance peaks depend on the electron number, as 
well as on the spatial symmetry. The results for large electron 
number and asymmetric potential are surprising: the shape of 
the peak-spacing distribution is Gaussian-like, the even/odd 
effect vanishes, and there is a substantial fraction of large spin 
ground states (S > 1). These effects imply a strong effective 
or residual electron-electron interaction. 

This is remarkable considering that conditions in our dots 
are not extreme: r s ~ 1.5 corresponds to a moderate bare 
interaction strength, and the dimensionless conductance is 
large, g ~ 4. In fact, previous work using the RPA-RMT 
approach suggested that we should obtain a strong even/odd 
effect Jl4] , |l6| ]. Conversely, to obtain the spin and peak- 
spacing distributions that we find here from the RPA-RMT 
model requires an effective exchange constant of J s ~ 0.6A, 
larger than the maximum value possible in RPA. The origin 
of this unexpectedly large residual interaction is not presently 
known, and we leave it for future investigation. 

We close with a caveat and a comment: First, this work 
is based on the 2D local spin density approximation whose 
validity, though already verified for small parabolic QD's 
[ |2^ , 30 1, is not well tested for large non-parabolic quantum 
dots as studied here. Our result highlights the need to go 
beyond RPA-RMT and perform a real Fermi liquid theory 
study. On the other hand, note the recent experimental work 
in Ref. ^3| which, though still preliminary because of insuf- 
ficient data, indicates a high probability of large-spin ground 
states. The surprisingly large effective interactions found here 
suggest that more experiments should be a high priority. 
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